Observation of scale invariance and universality in two-dimensional Bose gases 
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The collective behavior of a many-body system near a continuous phase transition is insensitive 
to the details of its microscopic physics[l|. Characteristic features near the phase transition are 
that the thermodynamic observables follow generalized scaling laws [J. The Berezinskii-Kosterlitz- 
Thouless (BKT) phase transition^, 0] in two-dimensional (2D) Bose gases presents a particularly 
interesting case because the marginal dimensionality and intrinsic scaling symmetry [J] result in a 
broad fluctuation regime which manifests itself in an extended range of universal scaling behavior. 
Studies on BKT transition in cold atoms have stimulated great interest in recent vears[5l4lOl]. clear 
demonstration of a critical behavior near the phase transition, however, has remained an elusive goal. 
Here we report the observation of a scale-invariant, universal behavior of 2D gases through in-situ 
density and density fluctuation measurements at different temperatures and interaction strengths. 
The extracted thermodynamic functions confirm a wide universal region near the BKT phase tran- 
sition, provide a sensitive test to the universality prediction by classical-field theory [n], U2] an d 
quantum Monte Carlo (MC) calculations [131], and point toward growing density-density correlations 
in the fluctuation region. Our assay raises new perspectives to explore further universal phenomena 
in the realm of classical and quantum critical physics. 
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In 2D Bose gases, critical behavior develops in the 
BKT transition regime, where an ordered phase with 
finite-ranged coherence competes with thermal fluctu- 
ations and induces a continuous phase transition from 
normal gas to superfluid with quasi-long range order [3]. 
In this fluctuation region, a universal and scale-invariant 
description for the system is expected through the power- 
law scaling of thermodynamic quantities with respect to 
the coupling strength and a characteristic length scale [L?! 
[hi ], e.g., thermal de Broglie wavelength (Fig. la). For 
weakly interacting gases at finite temperatures, in partic- 
ular, the scale invariance prevails over the normal, fluc- 
tuation, and superfluid regions because of the density- 
independent coupling constant [l5j and the symmetry of 
underlying HamiltonianQ. 

In this letter, we verify the scale invariance and uni- 
versality of interacting 2D Bose gases, and identify BKT 
critical points. We test scale invariance of in situ den- 
sity and density fluctuations of 133 Cs 2D gases at various 
temperatures. We study the universality near the BKT 
transition by tuning the atomic scattering length using 
a magnetic Feshbach resonance [16[ and observing a uni- 
versal scaling behavior of the equation of state and the 
quasi-condensate density. Finally, by comparing the lo- 
cal density fluctuations and the compressibility derived 
from the density profiles, we provide strong evidence of 
a growing density-density correlation in the fluctuation 
regime. 

We begin the experiment by loading a nearly pure 
133 Cs Bose condensate of N = 2 x 10 4 atoms into a 



'Current address: Department of Physics, The Pennsylvania State 
University, University Park, Pennsylvania 16802, USA 



single pancake-like optical potential with strong confine- 
ment in the vertical (z-) direction and weak confine- 
ment in the horizontal (r-) direction (l7l Il8j. The trap- 
ping potential, V(r, z) = mui 2 r 2 /2 + mu> z z 2 /2, has mean 
harmonic trapping frequencies u> r — 2ir x 10 Hz and 
lu z = 2ir x 1900 Hz. Here, r denotes the radial dis- 
tance to the trap center and m is the cesium atomic 
mass. In this trap, the gas reaches temperatures as low 
as T = 15 nK and moderate peak chemical potential 
Mo < ^bT. The ratio hio z /fio > ftu^/fe^T ~ 6 indicates 
that the sample is deeply in the 2D regime with < 1% 
population in the vertical excited states. Here, h = h/2ir, 
h is the Planck constant, and fcs is the Boltzmann con- 
stant. The 2D coupling constant is evaluated according 
to g — ^/Sna/lz (la ], where a is the atomic scattering 
length and l z = 200 nm is the vertical harmonic oscil- 
lator length. We control the scattering length a in the 
range of 2 ~ 10 nm, resulting in weak coupling strengths 
g = 0.05 ~ 0.26. Here, the density-dependent correction 
to g[lj|[l9| is expected to be small and negligible (< 2%). 

We obtain in situ density distributions of 2D gases 
by performing absorption imaging perpendicular to the 
horizontal plane with a commercial microscope objective 
and a CCD camera[18( (see Fig. lb for sample images). 
About 50 images are collected for each experiment condi- 
tion, and the average density n and the density variance 
5n 2 are evaluated pixel- wise (see Methods). We obtain 
the radial density n(r) and variance 5n 2 (r) profiles (Fig. 2 
insets) by accounting for the cloud anisotropy and per- 
forming azimuthal averaging[17j. 

We obtain the equation of state n(/j,T) from the av- 
eraged density profile by assigning a local chemical po- 
tential /j,(r) = /xo — V(r, 0) to each point according to 
the local density approximation. Both T and [1q can be 
determined from the low density wing where the sample 
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FIG. 1: Illustration of scale invariance and universality in 2D 
quantum gases, a, Scale invariance links any thermodynamic 
observable at different fi and T via a simple power-law scaling. 
In a 2D Bose gas with coupling constant j< 1, atomic density 
n measured at different temperatures (red lines) can be scaled 
through constant \ijT and n/T contours (dashed lines). Near 
the BKT phase transition boundary (green plane), systems 
with different g = (71,(72— (blue planes) scale universally, b, 
In situ density measurements of trapped 2D gases provide cru- 
cial information to test the hypotheses of scale invariance and 
universality. Sample images at different scattering lengths a 
are obtained from single shot. 



is assumed normal and the density profile can be fit to a 
mean-field formula n(/x, T) = — A^J ln[l — expfa/ksT — 
gn\j B /n)]^, where XdB = h/ \j2mnkBT is the thermal 
de Broglie wavelength. 

We confirm the scale invariance of a 2D gas by first in- 
troducing the dimensionless, scaled forms of density h 
n\ 2 dB (phase space density), fluctuation Sh 2 
and chemical potential jj, = /i/ksT, and showing that the 
equation of state and the fluctuation satisfy the following 
forms: 



5n 2 \j B , 



h = F(fi) 
5n 2 = G{p), 



(1) 

(2) 



where F and G are generic functions. This suggests both 
energy and length scales are set solely by the thermal 
energy and the de Broglie wavelength, respectively. An 
example at g = 0.26 (a = 10 nm) is shown in Fig. 2. 
Here we show that while the original density and fluc- 
tuation profiles are temperature dependent (see Fig. 2 
insets), all profiles collapse to a single curve in the scaled 
units. At negative chemical potential /2 < 0, the sys- 
tem is normal and can be described by a mean-field 
model (dashed lines). In the range of < fl < 0.3, 
the system enters the fluctuation regime and deviation 
from the mean- field calculation becomes evident. Cross- 
ing from normal gas to this regime, however, we do not 
observe a sharp transition feature in the equation of 
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FIG. 2: Scale invariance of density and its fluctuation, a, 
Scaled density (phase space density) h = n\^ B as a func- 
tion of the scaled chemical potential fi = fi/ksT measured 
at five different temperatures: T = 21 nK (black circles), 
37 nK (red squares), 42 nK (green triangles), 49 nK (blue di- 
amonds), and 60 nK (magenta stars), and coupling strength 
g = 0.26. Mean-field expectations for normal gas (dashed 
line) and superfiuid (solid line) are shown for comparison. In- 
set shows the radial density profiles before scaling, b, Scaled 
fluctuation 5n 2 = 8n 2 \\ B at different temperatures. Dashed 
line is the mean-field calculation based on the fluctuation- 
dissipation theorem [20||. Solid line is an empirical fit to the 
crossover feature from which the critical chemical potential 
jl c is determined. Inset shows the radial fluctuation profiles 
before scaling. The shaded area marks the fluctuation region 
< p, < p, c . Error bars show standard deviation of the mea- 
surement. 



state. At even higher jx > 0.3, the system becomes a 
superfiuid and the density closely follows a mean-field 
prediction[12] h = 2nji/g + \n(2ng/TT — 2/i). We no- 
tice that the mean-field theory in the superfiuid limit 
also cannot accurately describe the system in the fluctu- 
ation regime. Transition into the BKT superfiuid phase 
is most easily seen in the scaled fluctuation Sh 2 , which 
crosses over to a nearly constant value due to the sup- 
pression of fluctuation in the superfiuid regime 2l| . In 
the density profile n, a corresponding transition feature 
can be found when one computes the derivative dn/djl, 
i.e., the scaled compressibility k, as suggested by the 
fluctuation-dissipation theorem discussed in later para- 
graphs and Fig. 4. Finally, our measurement suggests 
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that the validity of scale invariance extends to all ther- 
mal, fluctuation and superfluid regimes, a special fea- 
ture for weakly-interacting 2D gas es H which underlies 
the analysis of a recent experiment [22j. 

We associate the crossover feature in the density fluc- 
tuations Sn 2 and the scaled compressibility k with the 
BKT transition pll [23j|. To estimate the location of the 
transition point, we apply an empirical fit to this feature 
and determine the critical chemical potential jl c and the 
critical phase space density n c (see full Methods). Re- 
sults at different g in the range of 0.05 to 0.26 are shown 
in Fig. 3c-d and compared to the theoretical prediction of 
n c = ln(£/g) and fi c = (g/n) ln^/^gl, where £ = 380 
and £ M = 13.2 are determined from a classical-field MC 
calculation (Tlj|. Our results show good agreement with 
the theory, apart from a potential systematic error from 
the choice of the fit function, which can account for a 
down shift of 10% in the fit values of jl c and h c . 

Further comparison between profiles at different inter- 
action strengths allows us to test the universality of 2D 
Bose gases. Sufficiently close to the BKT critical point 
with \p, — p, c \ < g, one expects the phase space density 
shows a universal behavior[l2l|. 

ii-n c = H^—^), (3) 
9 

where H is a generic function. Here, density and chem- 
ical potential are offset from the critical values h c and 
p, c , which remove the non- universal dependence on the 
microscopic details of the interaction [121, llil] . 

To test the universality hypothesis, we rescale jj, to p,/g 
and look for critical values n c and p, c such that the equa- 
tions of state at all values of g display a universal curve in 
the phase transition regime (see full Methods). Indeed, 
we find that all rescaled profiles can collapse to a single 
curve in the fluctuation region —1 < (fi — fi. c )/g < and 
remain overlapped in an extended range of \ jl — jji c \/g < 2 
(see Fig. 3a), which contrasts the very different equations 
of state n(fb) at various g shown in the inset of Fig. 3a. 
Our result closely follows the classical-field prediction fl2l| 
and quantum MC calculations [l3| assuming strictly 2D 
mean-field contribution, and the fitting parameters: crit- 
ical density fi c and chemical potential jl c show proper 
dependence on g and are in fair agreement with the the- 
ory prediction [ill (see Fig. 3c-d). We emphasize that 
critical values determined from the density fluctuations 
(see Fig. 3c-d) match well with those determined from 
the universal behavior, indicating that universality is a 
powerful tool to determine the critical point from a con- 
tinuous and smooth density profile. Similar agreement 
with the theory on the critical densities has also been re- 
ported based on different experiment techniques [rl IR [Toj . 

Further universal features near the phase transi- 
tion can be revealed in the growth of the quasi- 
condensate (QC) density n q = v n 2 — Sn 2 across the 
phase transition [ill Il2l. [25l| . QC is a measure of the non- 
thermal population in a degenerate Bose gas. A finite QC 
density does not necessarily imply superfluidity, but can 
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FIG. 3: Universal behavior near the BKT critical point, a, 
Rescaled density profiles h — h c measured at various coupling 
strengths, g — 0.05 (green triangles), 0.13 (blue diamonds), 
0.19 (red circles), and 0.26 (magenta squares) . Inset shows the 
original equations of state n(/x). b, scaled quasi-condensate 
density h q = \/n 2 — Sn 2 at different interaction strengths. In 
both plots, MC calculations from Ref. (open circles) and 
Ref. [13| ] (a, open squares for g = 0.07 and open triangles for 
g = 0.14; b, open squares) are plotted for comparison. The 
shaded area marks the superfluid regime and the solid line in 
b shows the superfluid phase space density calculation [l2T]. c 
and d, critical values fl c and h c determined from the following 
methods: universal scaling as shown in a (see full Methods, 
red squares), density fluctuation crossover (see text, black 
circles), and MC calculation from Ref. [ll| (solid line). Ex- 
periment values coincide at g — 0.05 identically, as a result of 
our analysis (see full Methods). Error bars show the standard 
deviation of the measurement. 

be responsible for a non-Gaussian distribution observed 
in momentum spaceQ. QC is predicted to be universal 
near the critical point following [12] 

n q = Q(^-^), (4) 
9 

where Q is a generic function and fi q = n q X 2 B . 

We employ both of our density and fluctuation mea- 
surements to evaluate n q at various g. Adopting fl c deter- 
mined from the universal behavior of the density profile, 
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FIG. 4: Fluctuation versus compressibility. Scaled compress- 
ibility R = F'(fi) and scaled density fluctuation Sn 2 = G(/i) 
are derived from measurements at two interaction strengths, 
g = 0.05 (squares) and g = 0.26 (circles), each containing two 
different temperatures between 20 and 40 nK (solid and open 
symbols, respectively). Diagonal line shows the expectation 
of G = F' in the normal gas region. Solid line shows sup- 
pressed fluctuation G = F'/(l + z) with z = 2. The color 
shaded areas mark the regimes of normal (left), fluctuation 
(middle), and superfluid (right). 



we immediately find that all measurements collapse to a 
single curve in the range of \p, — £t c \/g < 2 with appar- 
ent growth of QC density entering the fluctuation region 
(Fig. 3b). The generic function Q we determined is in 
goo d ag reement with the classical-field[l2j and quantum 
MC[l3j] calculations with no fitting parameters. Both 
our density and fluctuation measurements show univer- 
sal behaviors throughout the fluctuation region where a 
mean-field description fails and confirm universality in a 
2D Bose gas near the BKT phase transition (l2l. Il3| . 

The generic functions we described in the previous 
paragraphs offer new avenues to investigate the critical 
behavior of the 2D gas. Following the framework of scale 
invariance, we compare the dimensionless compressibil- 
ity k = dh/dfi = F'(p) and the fluctuation Sh 2 = G(j2) 
extracted from the measurements at g = 0.05 and 0.26 
(see Fig. 4). In the normal gas regime at low phase space 
density (G(jl), F'{jl) < 3), a simple equality G = F' is 
observed. This result is consistent with the fluctuation- 
dissipation theorem (FDT) for a classical grand canonical 
ensemble [20], which gives /c^T^ = SN 2 , where N is the 
particle number in a detection cell. In the fluctuation and 
the superfluid regimes at higher phase space density, our 
measurement shows that density fluctuations drop below 
the compressibility G < F' . 

Natural explanations for the observed deviation in- 
clude non- vanishing dynamic density susceptibility at low 
temperature[26j and the emergence of correlations in the 
fluctuation region (2 7] . While the former explanation is 



outside the scope of this article, we show that the corre- 
lation alone can explain our observation. Including cor- 
relation, the compressibility conforms top?! [28| 

«(r) = AjJ J (5n{r)6h(r + r'))d 2 r' (5) 
= Sh 2 (r)(l + z), (6) 



where (...) denotes 

l+w(r) J[g (2) (r,r+r')-l]d 2 



ensemble average and z — 
— 1 is the relative strength of 



l+n(r) J v [gW (r,T+r')-l]d 2 r' 

correlation to local fluctuation 5h 2 [27} . Here t/ 2 ) is the 
normalized second-order correlation function [29j and v 
denotes the effective area of the resolution limited spot. 
When the sample is uncorrelated, we have z = 0; non- 
zero z suggests finite correlations in the sample. In the 
fluctuation region shown in Fig. 4, observing a lower fluc- 
tuation than would be indicated by the compressibility, 
with z approaching 2, suggests that the correlation length 
approaches or even exceeds our imaging cell dimension 
y/v ~ 2 /im. This observation is in agreement with the 
expected growth of correlation when the system enters 
the fluctuation region. Similar length scales were also ob- 
served in the first-order coherence near the BKT phase 
transition using an interfcromctric method 8] and near 
the superfluid phase transition in three dimensions (30j] . 

In summary, based on in situ density measurements at 
different chemical potential, temperature, and scattering 
length, we have explored and confirmed the global scale 
invariance of a weakly-interacting 2D gas, as well as the 
universal behavior near the critical point. Our results 
provide detailed description of critical thermodynamics 
near the BKT transition and offer new prospects to inves- 
tigate other critical phenomena near classical or quantum 
phase transitions. In particular, we present experimen- 
tal evidence of the growing correlations in the fluctua- 
tion region through the application of the fluctuation- 
dissipation theorem. Further investigations into the cor- 
relations will provide new insights into the rich critical 
phenomena near the transition point, for instance, criti- 
cal opalescence and critical slowing. 



Methods Summary 

Preparation and detection of cesium 2D Bose gases are 
similar to those described in Ref. [l8| . We adjust the 
temperature of the sample by applying magnetic field 
pulses near a Feshbach resonance to excite the atoms. 
We then tune the scattering length to a designated value, 
followed by 800 ms wait time to ensure full thermalization 
of the sample. 

Absorption imaging is performed in situ using a strong 
resonant laser beam, saturating the sample to reduce the 
optical thickness. Atom-photon resonant cross-section 
and atomic density are independently calibrated. Aver- 
aged atom number N% and number fluctuation 5Nf at 
i-th CCD pixel are evaluated pixel-wise based on images 
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taken under identical experiment conditions. The pho- 
ton shot-noise, weakly depending on the sample's optical 
thickness, is calibrated and removed from the measured 
number variance. We correct the effect of finite imag- 
ing resolution on the remaining number variance using 



calibration from dilute thermal gas measurements. The 
density fluctuation 6nf is obtained from the recovered 
atom number variance using SnfX^B = /A, which 
replaces the dependence on the CCD pixel area A by a 
proper area scale \ 2 dB (details see full Methods). 
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Methods 

A. Calibration of the atomic surface density and 
the atom number fluctuation. 

Detection of caesium 2D Bose gases is detailed in ref. 
[IH and the atomic surface density n of the 2D gas is 
evaluated with similar schemes discussed in Ref. [3l[ , 
where the resonant cross-section cr is independently cal- 
ibrated using a thin 3D Bose condensate with similar op- 
tical thickness and the known atom number-to- Thomas- 
Fermi radius conversion. The calibrated value of <To can 
be compared to that determined from the atom shot- 
noise amplitude in dilute 2D thermal gases, where the 
noise is evaluated using binned CCD pixels to remove 
finite resolution effects. For dilute thermal gases, we ex- 
pect SN 2 — N, where iV is the mean atom number; we 
compare the fluctuation amplitude to the mean and ex- 
tract the value of <jq. Two results agree to within 10% 
and the residual non-linearity in the density calibration 
is negligible. 

We evaluate the atom number variance SN 2 pixel- wise 
based on images taken under identical experiment con- 
ditions. The photon shot- noise contribution 5N 2 , which 
weakly depends on the sample's optical thickness nao, is 
calibrated and removed from the atom number fluctua- 
tion using SN 2 = (5N 2 /2)[1 + (1+ J^? )2 e" a °], where 

SNq is the photon shot-noise without atoms and 7 is the 
ratio of the imaging beam intensity to the saturation in- 
tensity. Both 5Nq and 7 are experimentally calibrated. 
We then correct for the effect of finite resolution on the 
number fluctuation (l7| by comparing the atom number 
variance in a dilute thermal cloud to its mean atom num- 
ber, using SN 2 = N, and applying this calibration to all 
fluctuations measured at lower temperatures and higher 
densities. 



B. Density-density correlation in the fluctuation 
measurement. 

In the fluctuation measurement, we determine Sn 2 
from the pixel- wise atom number variance using the for- 
mula Sn 2 X 2 B = 6N 2 /A, which replaces the dependence 
on the pixel area A by a natural area scale X 2 B . This 
definition, however, does not fully eliminate the depen- 
dence on the imaging resolution spot size v ~ (2 /jm) . In 



particular, when the density-density correlation length £ 
approaches or exceeds the resolution, the measured fluc- 
tuation can depend on the fixed length scale s/v, which 
can complicate the scaling behavior. However, we do not 
see clear deviation of scale invariance and universality 
within our measurement uncertainties (Fig. 2b and 3b). 
We attribute this to the small variation of the non-scale 
invariant contribution within our limited range of sam- 
ple temperature. Further analysis on the correlations and 
fluctuations is in progress and the result will be published 
elsewhere. 

C. Determination of the BKT critical values from 
the fluctuation data. 

We use a hyperbolic function y(jX) — s(p, — p, c ) — 
yj s 2 (jl — jl c ) 2 + w 2 to empirically fit the crossover fea- 
ture of the density fluctuation near the transition region, 
assuming 8n 2 (jl) = De v ^\ where the critical chemical 
potential jl c , the fluctuation in the superfluid regime D, 
the slope of the exponential rise s, and the width of the 
transition region w are fitting parameters. The critical 
phase space density is then determined from the density 
profile as h c — h(fl c ). Other choices of fit functions give 
similar results, contributing only small systematics from 
the choice of different models. 



D. Obtaining the universal function H(x). 

We use the density profiles in the inset of Fig. 3a to 
look for critical values n c and p, c such that the equations 
of state at all values of g collapse to a single universal 
curve H(x) = h(ji) — h c , where x = (ft, — fic)/g is the 
rescaled chemical potential. To do this, we take the pro- 
file measured at g = 0.05 = g r as the reference, evaluate 
H r [x) = n(g r x + /2 c , r ) — n c ,r using the critical values 
n Cj r and £b Cjr determined from the fluctuation crossover 
feature, and smoothly interpolate the data to make a 
continuous reference curve H r (x) in the range of |x| < 1. 
Using this model, we perform minimum chi-square fits 
to the profiles measured at all other values of g accord- 
ing to h(fi) = h c + H r { ^~^° ), with only h c and jl c as 
free parameters. This procedure successfully collapses 
all density profiles (see Fig. 3a), and is independent of 
any theoretical model. The resulting critical values h c 
and jl c are plotted in Fig. 3c-d. 



